Thermoelectric power of MgB 2 - x Be a 
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We investigated thermoelectric power S(T) of MgEa-ssBe-e (x = 0, 0.2, 0.3, 0.4, and 0.6). S(T) 
decreases systematically with x, suggesting that the hole density increases. Our band calculation 
shows that the increase occurs in the cr-band. With the hole-doping, T c decreases. Implication of 
this phenomenon is discussed within the BCS framework. While the Mott formula explains only the 
linear part of S(T) at low temperature, incorporation of electron-phonon interaction enables us to 
explain S(T) over wide temperature range including the anomalous behavior at high temperature. 



I. INTRODUCTION 



Since the recent discovery of superconductivity in 
MgB2,u both theoretical and experimental efforts have 
been made to understand its structural and electrical 
properties. Most theoretical works suggest that coupling 
of the er-hole with B-plane phcoion is the key ingredi* 
ent of the superconductivity.cn 4 ] Boron isotope effectEl 
and other experimental data showed that the material 
is in the intcmAodiate or strong BCS coupling regiem 

(a = 0.7 ~ i)jhBtaE3 

Chemical substitution experiment has drawn much at- 
tention due to the possibility of obtaining higher super- 
conducting transition temperature T c . In general, by re- 
placing Mg- or B-sites with other elements, physical pa- 
rameters such as lattice constants and carrier density 
change. Therefore, study of T c shift in well-controlled 
substitution samples provides a chance to understand the 
superconductivity in detail. For the Mg-site substitu- 
tion, sevexalxompounds such as Mgi_ 2 .Al x B2 have been 
studied.trEjlLj For ,thp-|B-site, C- and Be-substitutions 
have been reported.ll3iI3 

It is important to know how these substitutions change 
the carrier density and how the change is related with 
T c shift. Thermoelectric power (TEP) measurement 
is a useful probe of the carrier density. For example, 
TEP experiment on Mgi_ x ALB2 showed that the Al- 
substitution dopes electrons.EJ Theoretical calculation on 
Mgi_ x Al x B2 suggested that the doping occurs largely in 
the cr-hole pocket.E3 In the present paper, we report our 
results of TEP measurement and theoretical band cal- 
culation on MgB2_ :c Be :E samples where B is substituted 
with Be. In an earlier paper, we showed that in this 
compound, the MgB2 phase is maintained up to x — 0.6. 
Also the lattice constants and T c changed systematically 
with Be-substitution.U We find that the TEP decreases 



with x, opposite to that in Mgi_ x Al ;E B2, which suggests 
that hole is doped. Also, changes in the band structure 
indicate that the hole doping occurs mostly in the cr- 
band. Thus, Be-substitution brings about hole-doping 
into the boron plane. We consider the consequence of 
this effect on the transition temperature and arue that 
the hole density change plays only a second role in the 
T c shift. 

Meanwhile, it is well known now that the temperature 
dependence of TEP is one of the unconventional features 
of MgB2i At low T, TEP is linear in T which is normal for 
most metals, but crosses over to a sublinear behavior as T 
increases. There-have been many interpretations of this 
phenomenon.talHJ We find that the crossover exhibits a 
systematic change with the Be-substitution. Furthur, we 
show that TEP in the whole T region is explained by a 
single model function in which the electron-phonon in- 
teraction is explicitely taken into account. 



II. EXPERIMENTAL 

Polycrystalline samples used in this experiment were 
synthesized by a powder metallurgical technique using a 
high pressure furnace. Starting materials are fine pow- 
ders (-325 mesh) of Mg (99.8%, Alfa Aesar), amorphous 
B (99.99%, Alfa Aesar), and Be (99.+%, Alfa Aesar). 
Stoichiometric amounts of powders are mixed and pel- 
letized. The pellets are placed in a tungsten vessel with 
a close-fitting cap, reacted two hour at 850 °C under 20 
atm. of high purity argon atmosphere. 

In MgB2_ x Bei;;, the MgB2 phase is maintained up to 
x = 0.6. Their structuraLand superconducting properties 
were reported elsewhere.EH We found that as a function 
of Be-substitution, lattice parameters show decreasing a- 
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and increasing c-values, and transition temperature T c 
decreases as summarized in Table I. 

For TEP measurements, bar-shaped samples (with di- 
mensions of ~ 4x0.5x0.1 mm 3 ) were mounted on two 
resistive heaters. Gold wires were used for thermoelec- 
tric potential leads. Chromel-constantan thermocouples 
were used for the temperature gradient measurement. 
Sample ends and thermocouple beads were glued to the 
heater blocks by Stycast epoxy. In our measurement, 
thermopower from the contact wires was carefully cali- 
brated. 



III. RESULT AND DISCUSSION 

Figure 1 shows thermoelectric power S(T) of 
MgB 2 _ x Be a; samples. For MgB 2 (x = 0), S(T) is linear in 
T at low temperature above T c . At higher T (> 150 K), it 
crosses over to a sub-linear behavior. These features to- 
gether with the large jump at-X^-i^ 1-4 /«V/K) are consis- 
tent with the earlier reports.HLj As boron is substituted 
with Be, the linear slope decreases. Also the crossover 
temperature is lowered. At x = 0.6, S(T) changes sign 
at a low temperature. 

Thermal diffusion of carriers gives rise to the linear-in- 
T behavior in S(T) at low temperatures. Kinetic trans- 
port theory shows that S(T) = {n 2 k 2 B T/3e)a (E)/a(E), 
where the dc-conductivity <j(E) and its derivative a (E) 
are calculated at the Fermi energy.Eil In a single parabeUs 
band system, it is approximated to the Mott formulaE^'EHl 



S(T) = X b T 



f7T 2 k%\ 



(1) 



where ks is the Boltzman constant, ± corresponds to 
the carrier sign, and Ep is the Fermi energy relative to 
the band maximum (or minimum). However, MgB2 is a 
multi band system where a- and 7r-bands coexist. Also 
the Fermi surface is not spherical. Thus, Eq. (1) can be 
applied only approximately here. The positive value of 
Xb in pristine MgB 2 suggests that the dominant carrier 
is hole. The decrease of Xj, with x suggests that the hole 
density increases. To obtain Xb, we fit the linear part of 
S(T) with Eq. (1) as shown by the dashed lines. Since 
the data are not extrapolated to zero, vertical shifts were 
needed in the fitrij 

To understand the behavior of Xb more quantitatively, 
we calculated band structures of MgB2- a: Be a ; as shown 
in Fig. 2. Here we used the local density approxima- 
tion (LDA) with the linearized augmented plane wave 
method. To account for the random substitution of Be, 
the virtual crystal approximation (VCA) was employed. 
Details of the calculation method were described by Mehl 
et al.Ej In current work, the experimental lattice param- 
eters in Table I were used. Four bands contribute to the 
Fermi surfaces: two cr-bands with B p XtV character give 
the 2D hole-type cylinders around the T-A line and two 
7r-bands form the 3D honeycomb tubular networks. The 



latter consists of one electron-type at the if -point and 
another hole- type at the if -point. Note that the most 
prominent change with x is the hole increase in the cr- 
bands. Change in the 7r-bands is relatively small. These 
results tell us that the Be-substitution dopes hole into 
the cr-bands. 

It is interesting to note that T c decreases with the hole- 
doping, similarly to the electron-doped Mgi_ a; AL c B2. 
Theoretical studies show that the superconductivity arise 
from the cr-band holes coupled with the B-plane phonons. 
T c = 39 K in MgB 2 can be produced from the McMillan 
formular using the coupling constartt_\ = 1.01 and the 
Coulomb pseudopotential fj,* = 0.13.QQ Here A is propor- 
tional to the density of state (DOS) of the cr-bands at 
the Fermi energy, Nh(Ef)- From our band calculation, 
we find that N h (E F ) =0.22 (eV _1 per cell) at x = and 
N h (E F ) =0.24 (eV" 1 per cell) at x = 0.6. (the cr-bands 
are highly 2D-like and the DOS increases only slightly 
with hole doping.) Then A will increase proportionally 
to become 1.09 at x = 0.6, if we assume the other pa- 
rameters determining A do not change, (see for example 
Eq. (3) of Ref. 3) This yields T c = 45 K, which is in 
sharp contrast with the observed T c decrease. This sug- 
gests that the other parameters change significantly with 
the substitution and their effects, overcome the Nh(Ef) 
contribution. In another paperjill we dealt with this is- 
sue and showed that the lattice constant change is the 
primary cause of the reduced superconductivity. 

Now we consider the unusual behavior of S(T), i.e., 
the deviation from the linear dependence at high T. In 
previous works, it-has been attributed to the minor cat 
rier contribution,E3 the thermally activated transport ,E3 
and to the phonon-drag effect EJ Here, we consider ef- 
fect of electron-phonon interaction on S(T). According to 
Kaiser, the interaction contributes to enhance the TEP 
through modi%ing the carrier mass and thus the ther- 
mal diffusionrJ Taking this effect into account, Eq. (1) 
is rewritten as 



S(T) = {1 + \\ s (T)]X b T, 



(2) 



where Xb is the slope in Eq. (1), A is the electron-phonon 
coupling constant, and A S (T) is a function which repre- 
sents the T-dependent thermopower enhancement: 



*k B T' 



(3) 



Here, the normalized Eliashberg function a F{u>) con- 
sists of the phonon density of states F(u>) and the cou- 
pling constant a. G s {fuo /ksT) is a function associated 
with thermal population of phonons. For MgB2,we 
culated A S (T) using a 2 F{uj) reported by Liu et all 
and fit the data with Eq. (2). Here we used A and Xb 
as fitting parameters. For x > 0, the Be-substitution 
into the B-plane will change, probably significantly, the 
phonon structure. Thus a and F(uj) will depend on x. 
Since they are not known, we took the values of MgB 2 in 
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calculating X s (T) . Thus A and X^ we estimate for x > 
samples are under large uncertainties. 

Figure 3(a) shows the fit for MgB2 (solid line). The 
bare diffusion part {X^T) and the enhancement part 
(AA S (T) x X&T) are represented with dashed and dash- 
dotted line, respectively. The latter has a broad maxi- 
mum at T ~ 215 K. Inset shows calculated behavior at 
higher temperature. Fig. 3(b) shows the fitting results 
for MgB2_ 2; Be a ;. Note that the fit is reasonable except 
the small deviation for x — 0.3. 

In Fig. 4, we summarize the linear slope X^ obtained 
from our analyses. The fitting results using Eq. (1) and 
Eq. (2) are shown with the filled-circles and the trian- 
gles, respectively. Note that Xf, from the modified model 
is smaller than that from the bare diffusion model. This 
is due to the enhancement effect contained in Eq. (2). 
Also, we estimated Xf, from the band calculation (the 
dash-dotted line). Here the Ep in Eq. (1) was taken 
from the cr-hole bands, assuming contributions from the 
other bands are negligible. EH The band calculation re- 
sult is closer to the modified model result at low doping 
region, x ^ 0.3, which supports the importance of the 
electron-phonon interaction effect. At x — 0.3, Xf, from 
the two fits exhibit a sudden drop. This drop may be 
related to the. observed structural change in the same 
compositionO The incomplete agreement between the 
fit result and the band calculation result may come from 
complex effects not included in this work such as the 
multi-band contributions and anisotropic transport. 

Regarding the electron-phonon coupling constant 
(EPC), we obtain A = 0.90 for MgB 2 . This isjajrood 
agreement with the earlier reports of 0.7 ~ l.OBu'u For 
x > 0, A increases to 0.98(a; = 0.2), 1.31(0.3), 1.34(0.4), 
1.47(0.6). This result is quite unusual because, as T c de- 
creases with x, A is expected to decrease. Recently, evi- 
dences show that MgB2 has two gaps. In this case, EPC 
from transport measurement (= X tr ) is different from the. 
EPC which determines the superconducting T c (— A SC ).Q 
Thus, the increase of TEP A (= Xtr) does not necessarily 
contradict with the T c decrease. One should also keep 
in mind that the increase of TEP A may be simply an 
errorneous effect that arise from the uncertainties in a 
and F(u) for x > we mentioned aboverj 

Now, let us consider the sign change in S(T) of x = 0.6. 
Sign change in TEP is widely obseriied in many alloyed, 
system*, for example, Ag-Au alloyJ^J YBa2Cu307_5,E£l 
NbNa;,E3, etc. In YBa 2 Cu307_(5, the change is observed 
as oxygen deficiency S increases. In NbN^, S(T) is com- 
posed of the diffusive (positive in sign) and phonon-drag 
terms (negative in sign). At high T, the former is domi- 
nant while the latter prevails, at low T. In the intermedi- 
ate T , sign change occurs E3 It is tempting to interpret 
our observation similarly: the negative S(T) may cor- 
respond to the phonon-drag effect. However, note that 
in the pristine MgB2, the phonon-drag feature is not ob- 
served. Furthur, the feature, if any, should be suppressed 
with Bc-substitution because the phonon-drag generally 
disappears as randomness is increased. Origin of the sign 



change is thus remains for future study. 



IV. CONCLUSION 

/From the TEP measurement and band structure cal- 
culation on MgB 2 _ a; Be a ; (x = 0, 0.2, 0.3, 0.4, and 0.6), 
we found that the hole density increases with x in the 
(T-bands . Thus, the Be-substitution dopes hole into the 
boron plane. The fact that T c shifts in the same direction 
(lowering) as the electron-doped case suggests that car- 
rier doping is not the primary route to control the transi- 
tion temperature in MgB2. This result is consistent with 
the 2D nature of the cr-bands. Furthur, we showed that 
the anomalous behavior of TEP at high temperature can 
be explained. by taking the electron-phonon interaction 
effect into account. 
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FIG. 1. Themoelectric power S(T) of MgBa-^Be^ {x = 0, 
0.2, 0.3, 0.4, and 0.6). Dashed lines represent linear fits to 
the data. 



FIG. 2. LDA virtual crystal band structures of MgB2- a; Be :r 
for x = (top), 0.3 (middle), and 0.6 (bottom). The experi- 
mental lattice parameters are used. The horizontal reference 
at denotes Ef- 



FIG. 3. (a) Thermopower data of MgB2 (open circles) and 
fit with the modified diffusion model (solid line). Dashed- 
and dash-dotted lines represent the diffusion and enhance- 
ment part, respectively. Inset shows the calculated behavior 
at higher temperature, (b) TEP data of MgB2-iBe I (open 
circles) and the modified diffusion fit (solid lines). 



FIG. 4. The linear slope of MgB2- a; Be 2 ,. • : determined 
from the Eq. (1). A : from the modified diffusive fit in Eq. 
(2). Theoretical results from the band calculations are shown 
with dash-dotted line. Solid lines are for eye-guide. 



TABLE I. Physical properties of MgB2- a; Be :E .|-Transition 
width AT C is determined from 10-90% transition.!!] 



X 


a (A) 


c(A) 


Tc (K) 


AT C (K) 





3.084 


3.522 


38.4 


1.2 


0.2 


3.078 


3.540 
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8.4 


1.5 
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